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Abstract 

Assuming the conventional Casimir setting with two thick parallel perfectly conducting plates 
of large extent with a homogeneous and isotropic medium between them, we discuss the physical 
meaning of the electromagnetic field energy Wdisp when the intervening medium is weakly dispersive 
but nondissipative. The presence of dispersion means that the energy density contains terms of 
the form d[ue(u)]/dco and d\uju(oj)}/ du. We find that, as Wdisp refers thermodynamically to a 
non-closed physical system, it is not to be identified with the internal thermodynamic energy 
U following from the free energy F, or the electromagnetic energy W, when the last-mentioned 
quantities are calculated without such dispersive derivatives. To arrive at this conclusion, we 
adopt a model in which the system is a capacitor, linked to an external self-inductance L such 
that stationary oscillations become possible. Therewith the model system becomes a non-closed 
one. As an introductory step, we review the meaning of the nondispersive energies, F, U, and W. 
As a final topic, we consider an anomaly connected with local surface divergences encountered 
in Casimir energy calculations for higher spacetime dimensions, D > 4, and discuss briefly its 
dispersive generalization. This kind of application is essentially a generalization of the treatment 
of Alnes et al. [J. Phys. A: Math. Theor. 40, F315 (2007)] to the case of a medium-filled cavity 
between two hyper planes. 
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I. INTRODUCTION 



For many years, after its discovery in 1948 [1| the Casimir effect was a theoretical curiosity, 
although it had evident applications to van der Waals forces and models of hadrons 
[3|. The Casimir formula for the quantum vacuum force between conducting plates was 
generalized to dielectrics by Lifshitz in 1956 jij, and the 1973 experiment of Sabisky and 
Anderson testing the Lifshitz prediction to a good accuracy, is well known. 

But the renaissance in studies of the Casimir effect began in 1997 with the work of 
Lamoreaux [6J. He measured the Casimir force between a conducting plate and a spherical 
lens, which, through the proximity force approximation 0, H, @] , agreed with expectations at 
something like the 5% level. (The accuracy of this measurement remains under some dispute, 
because various corrections, such as the effects of surface roughness, patch potentials, and 
finite conductivity, were not adequately taken into account.) In subsequent years, a variety 
of experiments were carried o ut, some of much greater accuracy and at considerably shorter 
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distances, Refs. [10|, [III, [12 
corrections [2~o| . 

One reason why the Casimir effect has attracted so much attention in recent years is the 
question of the temperature correction to the attractive force between real metal surfaces. 
At large distances (some micrometers) the relative thermal correction achieves several tens 
of percent, but at these distances the force itself becomes weak, and the experimental tech- 
nique is not sufficiently sensitive to give clear-cut results. At small distances, around 100 
nanometers, the measurements are claimed to be of high accuracy, about 1%, but at such 
distances the thermal correction is relatively small. On the theoretical side, the process 
of extracting the temperature dependence was carried out with the prescription given in 
Ref. [2l|. Because of inaccessibility of the effect to precise experiments, the issue was not 
reconsidered until the modern era, when Bostrom and Sernelius 22| recognized that this 
prescription could not be correct, and that necessarily the transverse electric reflection co- 
efficient at zero frequency must vanish for metals. This led to a reduction by a factor of two 
in the prediction for the slope of the linear high-temperature behavior (which would only 
be visible in experiments carried out at several microns), but it would predict a new linear 
temperature term at low temperatures, resulting in a 15% correction to the result found by 
Lamoreaux. Lamoreaux believes that his experiment could not be in error to this extent 
(23)]. Mostepanenko and collaborators have insisted that this behavior is inconsistent with 
thermodynamics (the Nernst heat theorem), becaouse it would predict that the free energy 
has a term linear in T at low temperature. Such a behavior is predicted by the so-called 
modified ideal metal model and also, as advocated by these authors, when the Drude model 
is applied to the case of a metal with perfect crystal lattice without impurities, in which 

m, SB], [26 



relaxation of conduction electrons is only due to scattering on thermal phonons 
Moreover, they assert that the precision Purdue experiments 18| . performed at T = 300 
K rule out the large thermal corrections predicted by the use of the Drude model with 
lattice imperfections taken into account 27j. The first Purdue experiment was performed 
at distances larger than 260 nm. More exact repetitions of that experiment j28| have been 
performed, down to 160 nm. We and others have responded that real metals do not exhibit 
this thermodynamic inconsistency, and that most probably the experiments are not so ac- 
curate as claimed @, H M. The situation is summarized in recent reviews BUI. In 
particular, the lack of a thermodynamic inconsistency has been conclusively demonstrated 
MS, by showing that the free energy for a Casimir system made from real metal plates 
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with impurities has a quadratic temperature dependence at low temperature. Further evi- 
dence for the validity of the notion of excluding the TE zero mode for metals comes from 
the recent work of Buenzli and Martin [35] , corroborating earlier work by these authors and 
others 36, 3?|], who show from a microscopic viewpoint that the high-temperature behavior 
of the Casimir force is half that of an ideal metal, a rather conclusive demonstration that 
the TE zero mode is not present. 1 

Our purpose with the present paper is not to study the temperature corrections in further 
detail. The brief survey above indicates that there is a need of reconsidering the underpin- 
nings of the Casimir theory in some detail. As an attempt to do this, we will consider the 
Casimir problem from an unconventional angle, emphasizing the role of dispersive media. 
Thus, in Sec. Ill II we will show how the Casimir energy Wdisp for a dispersive nondissipative 
medium, reflecting a non-closed physical system, is not to be identified with the internal 
thermodynamic energy U, or the electromagnetic energy W, when U and W are calculated 
as though dispersion were not present. In this regard a capacitor model of the system proves 
to be quite illuminating. Finally, in Sec. [IV] we examine another aspect of phenomenological 
electrodynamics, namely its generalization to higher spacetime dimensions, D > 4, both 
because the topic has some relationship to the dispersive theory discussed in Sec. 1111] and 
also because the higher-dimensional electrodynamical theory is a topic of general current 
interest . 



II. FREE ENERGY F, INTERNAL ENERGY U, AND ELECTROMAGNETIC 
ENERGY W 

In order to fix the notation, and for reference purposes, we give in this section a brief 
survey of how the various energy concepts occur in Casimir theory. Assume the usual 
configuration, in which there are two thick infinitely large plates situated at z = and 
z — a, with a homogeneous and isotropic medium in between. We take this intervening 
medium to have permittivity e and permeability \i. In this section we take these material 
parameters to be constants. For simplicity we assume that the medium to the left (z < 0), 
as well as that to the right (z > a) are ideal (the so-called IM model), so that the TE and 
TM reflection coefficients t~te and ttm satisfy r^ E = r\ M = 1 for all Matsubara frequencies 
m, including m = (the breakdown of this assumption for the t"te coefficient at m = is 
the crux of the temperature controversy for real metals). Let n = y^ejl be the refractive 
(temperature independent) index of the intervening medium, (3 = 1/T, ( m — 2irm/j3, and 
k, 2 = k\ + n 2 ( 2 . The free energy F per unit surface area can now be written 

F = —;J2' koIk\ii (I - e- 2Ka ) , (2.1) 

for arbitrary T. The prime on the summation sign means that the m = term is counted 
with half- weight. 

The internal energy per unit area U is now constructed from the thermodynamical formula 



1 It could here be added, as a contrast, that Intravaia and Henkel have recently claimed that for metals 



with perfect crystal lattices the Lifshitz theory leads to violation of Nernst's theorem 38 1 
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From Eq. (12. ip it is apparent that f3 appears only in the lower limit of the integral in the 
expression for (3F. Since d( m /d/3 = —27im//3 2 , we get 

A 2 00 

^ = ^5> 2ln (l-e-« m ), (2.3) 

P m=0 

where 

Anna m . /N 

a = — — = AirnaT. (2.4) 

P 

The m = term does not contribute. One way of processing the expression (12.31) is to 
expand the logarithm, 

00 

m 2 ln(l -e~ Qm ) = -^-mV Qfcm , (2.5) 

k=l 

and then sum over m, whereby we get 

U = f 1 CQth 2 (27T " maT) . (2.6) 

^ m sinh 2 (27rnmaT) 

When n = 1, this agrees with Eq. (18) of Ref. [40j. (Cf. also the discussion of energy 
and free energy in Ref. j4lj.) The expansion (12.61) is most convenient at high temperatures, 
aT 3> 1. By including only the m — 1 term, one gets 

U = -A7rn 2 T 3 e- 4nnaT , aT > 1. (2.7) 

It is apparent that U — > when T — > 00. This is as we should expect physically: The Casimir 
energy measures the change in energy induced by the boundaries, and these constraints 
decrease in importance when the classical thermal energy becomes high. 

To get a convenient expression at low T one may perform a Poisson resummation, along 
the same lines as discussed in Ref. Define the quantity b(m), 

b(m) =m 2 ln(l-e-° H ) , (2.8) 

along with its Fourier transform c(q), 

c(q) = — [ b(x)e' iqx dx. (2.9) 
Then, according to the Poisson formula, 

b(m) = 2rr c ( 2rrm ) = 2 / x 2 ln(l- e~ ax ) dx 

m=~oo m=—oo ^ 

00 «QO 

+ 4^ / x 2 cos(27rmx)\n(l - e' ax ) dx. (2.10) 

m=l ^° 

The various terms can be evaluated analytically. The following formulas are here useful 
(Ref. 0, sec. 3.951), assuming b > 0, 



00 x 2m smbx , , xm <9 2m r 

00 x 2m+1 cos bx , , xm <9 2m+1 
rfx = (-1) 



, , 1 

— COtil7T0 

2 26 



^ , , 1 

7— — — COth 7TO 

db 2m+1 2 2b 



(2.11a) 
(2.11b) 
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The calculation gives, for arbitrary T, 



U = 27rn 2 T 3 



7T 



1440(naT) : 



+ 



naT 



oo 1 

Y — 



tt° * — ' rrr 

771=1 



nm irm 
coth ■ 



2naT 



2naT 



+ 



I _7rm_ V 
V 2naT / 



7rm , 7rm 
1 + - — - coth 



2naT 



2naT. 



It is of interest to consider the limit of low dimensionless temperatures, 

nal \ 



(2.12) 
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7T" 



720na 3 



1-720 



7T 



-1 C(3) + 48(naT) 4 



aT < 1. 



(2.13) 



Again, this agrees with the low-temperature expression obtained earlier, for instance in 
Ref. [43(, when n = 1. It is to be noted that U, as well as the corresponding low-temperature 
expression for F, 



7T 



72(W 



1 + 360 



naT 



7T 



C(3) - {2naTf 



(2.14) 



contain a term that is independent of a, which means that this term does not contribute to 
the force between the plates. 

The third kind of energy that we shall consider is the electromagnetic energy W, still 
taken per unit surface area. As above, we take the medium to be nondispersive. We start 
from the energy density, 



™ = \e{El + E])+ l -^Hl + Hl), 



(2.15) 



so that, per unit area, W = wa. Quantum mechanically, the product E^(r) is to be replaced 
by the expectation value {E z (t)E z (t')) in the limit when r' — > r. Similarly for the other 
components. We assume first that T = 0. According to the fluctuation-dissipation theorem 
in Fourier space we have 



i(Ei(r)E k (r')) w = lmT ik (r,r';uj), 

i{Hi(r)H k (r')) w = -r— ^curl^ curl^ ImTj^r, r'; to) 
H uj 



(2.16a) 
(2.16b) 



where curljfc = eijkdj, €ijk being the Levi-Civita symbol. Further, T is the Green's function 
as defined by Schwinger et al. [2l|], in terms of a polarization source P, 



with 



E(x) 



T(x,x') 



d i x'Y(x,x')-¥(x') 



duj 
2^ 



e~ luJT T(r, r'; u) 



and r — t — if . Introducing a transverse Fourier transform, 



(2vr) ; 



e ik x .(r-r') Z ';k ± ,u) 



(2.17) 
(2.18) 

(2.19) 



we can write 



K 1 

9xx = ^ cosn K ( z ~ z ')i (2.20a) 

& = —\ cosh K ( z - z ')i ( 2 - 20b ) 



where 



k 2 1 

g? = cosh k(z-z'), (2.20c) 



d = e 2Ka - 1, = fci - raV. (2.21) 



(Details are given in Ref. |44j.) (Note that the notation is slightly different than that given 
in Ref. gj.) 

Defining the Fourier components (..) w k of the energy density according to 

we first obtain for the electric part, letting z' — > 2, 

= j§(fl£ + ^ + «£) = (2-23) 

Then defining the magnetic counterpart gf k to the electric part gf k according to 

1 

we obtain by an analogous calculation, in the limit when z' — > z, 



9?k = — curliz curl' fcm gf m , (2.24) 



\»{H 2 U = + fl£ + ifi) = (2-25) 

The electric and magnetic contributions to the energy are equal, as we would expect. Adding 
the expressions (12.231) and (I2.25j) and multiplying with a we obtain, at zero temperature, 

where a frequency rotation u — > i( has been performed. This expression can be further 
processed by introducing new coordinates X = k± = kcosO, Y = n( = K,sm8, with 
k = \fk\ + n 2 ( 2 . We get 

1 r°° Z 3 dz TV 2 
~ ~487r 2 na 3 J e z - 1 ~ ~720na 3 ' 

in accordance with Eqs. (12.131) and (12. 14[) . 
At arbitrary temperature we get 



W = -8nn 2 aT 3 £ ™ 2 I (2.28) 



m=l 
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with k = \fk\ + (2irnmT) 2 . Alternatively, we may write 

W = -4vrn 2 T 3 ^ m 2 / (2.29) 

™_i Jam ^ ^ 



m=l 



where a = AnnaT as before. 

At high temperature, aT ^> 1, it is easy to check that W agrees with U calculated 
previously. We approximate the integral in Eq. (12.291) by J^ n e~ z dz = e~ am , and so get 

00 

W = -Ann 2 T 3 m 2 e" am -> -47m 2 T 3 e - 4 ™ aT (2.30) 

m=l 

when m — 1, in agreement with Eq. (12.71) . 

We shall not delve further into a detailed study of the equality between W and [/ in the 
case of arbitrary T. The equality should be clear on physical grounds, since we are dealing 
with a closed thermodynamical system. 

After having given this survey, we have the necessary reference background for studying 
the dispersive regime. 



III. ON THE DISPERSIVE CASE, NEGLECTING DISSIPATION 

As mentioned, our main focus will be on the dispersive case. Assume first that the 
medium in the region < z < a is both electrically and magnetically frequency dispersive, 
e = e(cj),n = fi(u). The walls are taken to be perfectly conducting, as before. The total 
energy density w<n S p is known to be |46|, 147 



Wdisp 



duj 
2^ 



d 2 k A 

W) 



d(euj) 
du 



<£ 2 U + 



d(fiuj) 
du 



(3.1) 



We can write this as a sum of two parts wj and wu, where wj is the same expression as in 
Eq. (I2.22p with e — > e(lu), \i n{w), and where 



w u 



dui 

£ 

2vr 



-uj 



d 2 k A 

W) 



de 
du 



+ ^(H 2 } 
duo 



(3.2) 



Correspondingly, for the surface densities, Wdisp = Wi + Wu. 

The first property to be noted in connection with Eq. (13.11) is that it is derived under 
the assumption of negligible dissipation. Some dissipation is always present — this being a 
consequence of Kramers-Kronig's relations — but it is a legitimate approximation to neglect 
it except in the neighborhood of eigenfrequencies in the cavity. This assumption means that 
the relaxation frequency in the dispersion relation can be set equal to zero, and we may 
adopt the usual dispersion relation for a dielectric, for simplicity taking it hereafter to be 
nonmagnetic, 

e ^ = 1+ r^s- " = L (3 ' 3) 

In the case of a general dissipative medium, neither the energy nor the stress tensor are 
derivable in terms of permittivity/permeability alone, and therefore cannot be given in a 
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general form using macroscopic methods. (This point is discussed in detail by Ginzburg 

Second, it is clear that the expression (13. ip is not intimately related to the Casimir effect 
as such. It is more natural to consider the problem as belonging to classical electrodynamics, 
namely a system of two conducting plates between which there are stationary electromag- 
netic oscillations. The expression (13.11) is actually obtained from the energy balance equation 

V(ExH) + ED + HB = 0. (3.4) 



(See, for example, Eq. (7.5) in Ref. |47|.) In order to accumulate electromagnetic energy, 
one has to consider oscillations that are not purely monochromatic, but distributed within 
a band of frequencies around each eigenfrequency. In this way external agencies, outside of 
the plates, are called for. It is natural here to regard the system to be a capacitor, linked 
to an external appropriately adjusted self-inductance L such that stationary oscillations 
become possible (external resistances are forbidden since we omit dissipation). That means, 
the plates with the intervening medium is thermodynamically a non-closed system. From 
this we can draw the important conclusion that the full dispersive energy M^disp is not to 
be identified with the thermodynamical energy W = U calculated earlier. The laws of 
thermodynamics are applicable to closed systems only. 

The mentioned model of a classical electromagnetic non-dissipative circuit is studied in 
Ref. [46}. It is instructive to consider the salient features of the argument also here: 

Let the charges Q be supplied and withdrawn from the plates with frequency to. The 
self-inductance of the circuit is L, as mentioned, and the electromotive force we call S. The 
potential across the plates is determined from the equation 

<p = £-LJ, (3.5) 

where J = Q. The frequency of the circuit is 

uj = l/y/LC(u), (3.6) 

where the capacitance C(oS) of the capacitor is determined by <ft — Q/C{uj). By considering 
almost monochromatic quantities [the same kind of argument that led to Eq. (13.11) ]. we get, 
when taking the average over a period, 

ei^AUr^mA. (3.7, 



dt [2 2 duo 

The expression between brackets is the circuit energy. From J = —iuoQ and Eq. (13. 6p we 
get \LJ 2 = \C cj) 2 and so the circuit energy may be written 

W*. = (3.8) 
2uj du 

This expression, because of the derivative with respect to a;, is seen to be related to Eq. (13.11) . 

Now consider a small adiabatic displacement of the capacitor plates. As W c \ tc /uj is an 
adiabatic invariant, 

6W ciTC = W^ c 8u/u. (3.9) 

By means of Eq. (13. 6p , 

5u 15C . 



s 



The change in C consists of two parts, 



dC 

8C=(5C) st + —5id, (3.11) 
did 

where the first term is the static part and the second term depends on the frequency change. 
From Eqs. flBTTUj) and flBTTTj) . 

, C<t = _i_^ fo . (3 . 12 ) 
id 2 duo 

When Eq. (13.81) is substituted in Eq. (I3.9P and (13.121) is used, dC/cLu disappears, and we get 

6W cilc = ~W{SC) Bt = - l -^(5C) st . (3.13) 

This is the same expression as one obtains by taking the variation of the average of the 
energy Q 2 /2C of a thermally insulated capacitor. It means that when dispersion is present, 
the electromagnetic stress tensor contains no derivatives with respect to the frequency. The 
argument is general, and is not critically dependent on our choice of a capacitor model. 
When applied to our case, we can thus conclude as follows: 

1. The dispersive energy Wdisp whose density is given in Eq. (13.11) refers thermodynam- 
ically to a non-closed system, and is therefore not to be identified with the internal 
energy U calculated in Sec. [Til starting from the free energy F, or the electromag- 
netic energy W, in the nondispersive case. As was demonstrated, when e and \i are 
constants, W = U . We are still to use the same expressions for W and U when the 
permittivity and permeability depend on frequency. 

2. As for the electromagnetic stress tensor, the derivatives with respect to id are not to 
be included. That is, the electromagnetic force can be calculated from Eq. (I2.22p with 
e — > e(cj), ji — > h(oj). 



It may finally be noted that by inserting the simple form (13.31) for e{id) for a dielectric, we 
obtain for the dispersive correction Wu = awn a divergent expression, 

— ^3 — 7o ^mwJ# u (3 - 14) 

cf. Eq. fl£2p . 

Another way to to see that the dispersive medium should be treated without the frequency 
derivative of the permittivity is to recognize that the Casimir energy may be derived by a 
variation expression 

SE i f did d 2 kj_ . N F , . . 

dzde{z)g kk (z,z,k ± ,id), (3.15) 



A 2 J 2tt(2tt 



which is Eq. (2.26) of Ref. [2l[. This starting point is equivalent to the variational argument 
recounted in this section. 
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IV. DISCUSSION ON AN ANOMALY IN THE CASIMIR ENERGY FOR 
HIGHER DIMENSIONS 



The electromagnetic theory in a continuous medium has in general a rich structure. 
Most notable is the fact that (within the commonly accepted Minkowski theory) the spatial 
photon momentum is equal to k = rau;k, implying that the photon four-momentum becomes 
spacelike, k^k^ = {n 2 — 1)uj 2 > (we make use of the metric g^ v = diag(— 1, 1, 1, 1), and 
assume n constant to begin with). Accordingly, there are inertial systems in which the 
photon energy becomes negative. A striking example of this sort is provided by the Cerenkov 
effect: In the inertial system where the emitting particle is initially at rest, the recoil kinetic 
energy of the particle is necessarily positive. Thus, in order to preserve energy conservation, 
the energy of the emitted photon has to be negative. Another example of a related sort 
is provided by the so-called anomalous Doppler effect, occurring when a quantum particle 
detector moves superluminally in the medium. Thus Ginzburg and Frolov [49|] studied 
such kinds of particle detectors and showed how the excitation of a detector uniformly 
accelerated in a vacuum with the associated emission of radiation is actually similar to the 
radiation occurring in the region of the anomalous Doppler effect when the detector is moving 
superluminally with constant velocity in the medium. See also the discussion in Ginzburg's 
book [il]. Situations of these kinds were discussed also by Brevik and Kolbenstvedt, in the 
case of constant velocity [Ho[ and for constant acceleration [5l| . 

We shall round off our paper not by considering the above-mentioned effects any further, 
but instead another effect that has also a bearing on medium electrodynamics, namely the 
anomaly that turns up in the case of higher spacetime dimensions, D > 4. The anomaly 
reflects the breaking of conformal symmetry. We do this because the topic has some rela- 
tionship to that considered in Sec. IHH and also because it has attracted interest recently 
in the case of a vacuum field. A generalization to the medium case thus appears natural. 
Higher dimensions, in the context of Casimir theory, were considered long ago by Ambj0rn 
and Wolfram 52 L but anomalies of the type considered below were not studied until recently 
by Alnes et al. 3f|; cf. also Refs. 53, 5J]. 

Let us assume, then, that there are two parallel hyperplanes with separation a, the region 
< z < a being filled with an isotropic medium of refractive index n = ^/ejl. The walls are 
assumed perfectly conducting, as before. The anomaly we wish to consider is present also 
in the case of zero temperature, so we shall assume T = in the following. 

The appropriate electromagnetic energy-momentum tensor is the Minkowski expression, 
called S^, 

Sffu = F ^H y a - -g^ u F af3 H al3 ; (4.1) 
cf., for instance, Refs. 55, 56|, 57]. Here F^ u = d^A u — d u A^ with /i, v = 0, 1, 2, ...D — 1 = d is 



the field tensor, whereas F ok = Ef. with k = 1, 2, ...d are the components of the <i-dimensional 
electric field vector E. The magnetic induction (B in the three-dimensional case) does not 
constitute a vector in the higher- dimensional case, but is given by the d(d — 1)/2 components 
of the antisymmetric spatial tensor F^. Analogously, the second D-dimensional tensor 
H^y occurring in Eq. (14.11) is given by the vector components H ok = D^, D being the d- 
dimensional induction vector, and by the d(d — l)/2 components of the spatial magnetic field 
tensor (H in the three-dimensional case). In analogy with three-dimensional theory, we 
assume constitutive relations in the form H oh = eF oh and F^ = [lH^ also when D > 4. 
Turning now to physical quantities, it is convenient to start with the surface pressure P 



10 



on the hyperplane z = 0. We observe that the usual expression for P (cf. Eq. ( 12. ip ) can 
easily be generalized to the case of d = D — 1 spatial dimensions. Taking into account that 
there are (D — 2) physical degrees of freedom in the field in the cavity, we have 



F = (D -2) 



2tt 



d d ~ l k 



(27T) 



-2na 



where 



2/-2 



k± — + ky + 



4- k 2 



(4.2) 



(4.3) 



The volume element in momentum space is d d x k± = 2 dk±, where the solid angle 

is determined by Qd-i — 2-K d ^ 2 [T(d/2)]~ 1 . The pressure P = —dF/da can now be written 



2(D 



(27T) 



&d-2 



d( 



tzfcj i dk i 



D 2Ka 



(4.4) 



The double integral over ( and can be further processed by introducing polar coordinates, 
2l| . Again, we introduce X = kcosO = k±, Y = nsin6 = n(, satisfying X 2 + Y 2 



K 2 . 



The area element in the XY plane is ndndQ = ndk±d(. The integral therewith becomes 



d( 



nk d r 2 dk i 



2ko, 



- 1 



n/2 



COS 



d-2 



ede 



K d dn 

2na 



(4.5) 



We now make use of known integral formulas and insert the expression for fld-2, to get 

(!>-2)(£>-l) T(D/2)((D) 



P 



n 



(4%) D / 2 a D 



(4.6) 



It ought to be emphasized that this expression was obtained without any regularization 
procedure. The presence of the medium is seen here to turn up through the factor n in the 
denominator. If n = 1, including the case of a vacuum as well as the case of a "relativistic" 
medium satisfying e = the expression reduces to that derived earlier [521 ] . This result 
parallels that obtained in the T = parts of the energy, cf. Eqs. (I2.13P and (12.141) . It is 
nearly identical to the result found in Ref. [45j for the scalar case in D dimensions [Eq. (2.35) 
there], differing only in the factor (D — 2)/n. 

The electromagnetic field energy density w in the cavity is a more delicate quantity. The 
natural way to calculate w is via the energy-momentum tensor. This procedure - carried out 
by Alnes et al. in the case of a vacuum cavity 139, l54j - led in the case of metallic boundary 
conditions to the result 



w 



(D — 2)r(Z?/2) 



where 



{47r) D / 2 a D 
Id (—) = s// ( D. — ) -h ' s u [ D. I - 



Wi + w 2 , 



(4.7) 



(4.8) 



(h being the Hurwitz zeta function. Note that the first term yields the pressure ( 14.61) . 

d 



— aw 1 
da 



(4.9) 
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so that the second term in the energy density, W2, which diverges like z~ D close to the surface 
when D > 4, does not contribute to the force between the plates. It can explicitly seen that 
written in physical variables this term is independent of the separation between the plates 
and hence does not contribute to the force. This anomaly can actually be seen to manifest 
itself in another way if we go back to the expression (14.11) for the energy-momentum tensor: 
its trace Sjf^ is nonvanishing when D > 4. Physically, as emphasized in Ref. 3{|, the 



divergent self energy of a single surface is related to the lack of conformal invariance of the 
electromagnetic Lagrangian for D > 4. All of this is exactly as seen in Ref. (45[, Chap. 11, 
for the scalar field. 

It turns out that the anomaly can be regularized away by subtracting off the self energy 
for both plates. Then, the second term in Eq. (14.71) is absent, and only the first, finite, terms 
in w remains. 

As mentioned, these calculations of w were made for the case of a vacuum cavity. The 



result was found via a combination of dimensional and zeta function regularizations [58 
The result could be recalculated for a medium cavity, but such a detailed calculation is 
hardly justified in view of the simple occurrence of n in the expression (14. 6p . In fact, since 
in physical units, w = TicjaP times a function of D, it is clear (for example, Ref. (47| . 
Eq. (36.12)) all we have to do to insert a uniform medium between the plates is replace c 
by c/n; this shows that the same factor n will appear in the denominator of the expression 
for w as it did in W or P. Thus, after regularization, we obtain the relationship 

P=(D-l)w 1 , (4.10) 

which is the same connection as for a vacuum. 

Finally, we consider an alternative method for obtaining the energy W that avoids the 
field theoretical approach above, and instead starts by considering the individual photon 
momenta directly. The photon momentum in the medium is \fk\ + 7r 2 m 2 /a 2 , and the 
photon energy is obtained by dividing this expression by n, assuming that n = constant. 
Thus we have, still at T = 0, 

7(2-1 



1 

W = ~ £ / \A* + - 2 -Va 2 - (4.11) 

In order to extract a finite expression we have to regularize in some way, for instance by 
using an exponential cutoff. The important point in our context is however that the integral 
in Eq. (j4.11j) . and the sum, are just the same as in a vacuum field. Thus the influence of 
the medium turns up only in the prefactor 1/n, in accordance with what was found above. 

Can this theory be generalized to the dispersive case? Such a performance is not quite 
straightforward, in view of the complicated form (13.11) for the dispersive energy density. 
Some insight can however be obtained from the following argument. As noted in Sec. IHfl 
the electromagnetic stress tensor does not contain derivatives with respect to the frequency. 
Thus, the momentum flux density has the same form as in a nondispersive medium. We 
may assume therefore that the photon wave vector is equal to n(u)u: 

Jtf_ + n 2 m 2 /a 2 = n(u)uj. (4.12) 

When the wave vector is given, this equation can be solved (numerically) for n{uj) and u. 
Inverting, we find n as a function of \fk\-\- ir 2 m 2 /a 2 . Accordingly, we can write the energy 
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as 



w = y f d"- 1 k± Vkj + vrW/a 2 

^ I (27r) d - 1 n(y/k 2 ± + vr 2 m 2 /a 2 )' ' 



m=l 



For very high wave numbers, n — > 1, and the integral reduces for these frequencies to its 
vacuum counterpart. One should bear in mind that the expression (14.131) holds only in an 
approximate sense, as we have ignored the accumulation of energy during the slow building 
up of the electromagnetic field. 



V. CONCLUSIONS 

Some care ought to be taken when dealing with dispersive and dissipative media. In the 
general case of arbitrary dispersion (which implies necessarily dissipation also), the elec- 
tromagnetic energy cannot be rationally defined as a thermodynamic quantity at all. If 
dispersion is weak, making it possible to ignore the accompanying dissipation, it is mean- 
ingful to define the electric energy such that it contains terms of the type d[uje(u)]/du, and 
similarly for the magnetic field. In such a case, as we have seen, one should distinguish 
the electromagnetic energy W<a S p from the thermodynamic energies calculated for nondis- 
persive media and used without the derivatives on e and fi in the dispersive case where 
n = n(u). The main reason for the difference is that in the dispersive case we are dealing 
with a non-closed physical system. 

The electro dynamic theory of media, especially when dispersion is included, has a rich 
structure. As an example of this, we showed in Sec. [IV] the anomaly turning up when 
D > 4, which is especially interesting when dispersion is present. This shows the interplay 
between local surface energy divergences and the breaking of conformal symmetry. The 
clarity brought to bear by the above analysis will now allow us to understand more fully the 
dispersive case, and to some extent also the questions connected with temperature problems. 
Moreover, we hope to have contributed to the understanding of surface energies and their 
significance. 

Finally, we make the following comment. Our electromagnetic formalism in this paper 
has been the conventional one, whereby the basis for calculating stresses on matter is the 
Abraham-Minkowski stress tensor (cf., for instance, Ref. 56J). Now, in a recent paper Raabe 
and Welsch [59] have developed a somewhat unconventional theory for electromagnetic fields 
in a medium based upon the Lorentz force, from which they derive a stress tensor different 
from the Abraham-Minkowski form. The Casimir effect was chosen by these authors as the 
physical phenomenon to which they applied their proposed theory. We merely mention this 
novel formulation here; it would lead us too far from our main purpose to make a detailed 
scrutiny of this rather complicated formulation. A comment on some consequences of the 



altered stress tensor in practical applications is under preparation [60 
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